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The traditional calculation of dilcpton yields from bremsstrahlung relies on the as- 
sumption that electromagnetic and strong processes factorize. We argue that this 
assumption, embodied by the soft photon approximation, cannot hold true for invari- 
ant mass spectra on very general grounds. Deriving a formula for the dilepton cross 
section for pion-pion scattering that does not rely on such factorization, we formulate 
the problem exactly in terms of three-particle phase space invariants. Using a simple 
\Q . one boson exchange model for comparison, we find that dilepton cross sections and 

' yields are generally overestimated by the soft photon approximation by factors 2-8. 

In extreme cases, overestimation up to a factor 40 is possible. 

a: 
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1 Introduction 

Interest in the use of dileptons as a probe of the hot and early phases of heavy ion collisions 
is fed by the desire of finding new physics such as the vaunted quark gluon plasma and 
generally probing the behavior of nuclear matter under extreme conditions Q . Recent data 
by HELIOS and NA38 and discrepancies found by CERES § have provided impetus to 
hopes that new physics is finally in sight. Realising such hopes, however, requires detailed 
understanding of background processes: the contribution from each must be calculated 
quantitatively. 

While quantitative calculations are challenging already on a technical level, it is a much 
harder problem to identify untested assumptions that enter such calculations and to quan- 
tify their effects. We here aim to show, by example of bremsstrahlung from pion-pion 
scattering, that the "soft photon approximation" (SPA) represents such an untested as- 
sumption Q. 

1 To be published in International Conference on Structure of Vacuum and Elementary Matter, Wilder- 
ness, South Africa, March 10-16, 1996 (World Scientific). 
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2 Why the SPA must fail 



Conventional dilepton yields at very low invariant masses M < 500 MeV are dominated 
by Dalitz decays and bremsstrahlung. Calculations of such bremsstrahlung yields rely 
heavily on the soft photon approximation because it is simple to use. This simplicity is 
achieved mainly through the assumption that the electromagnetic and strong processes 
f actor ize. 

In order for this assumption to be valid, two conditions || must be met: First, the photon 
energy go must be much smaller than the energy Ei of any one of the hadrons participating 
in the scattering, qo/E{ <C 1; secondly, the hadronic and electromagnetic scales must be 
sufficiently different to permit separate treatment. This translates into the condition 
qo <C mylPil/Ei, where my is the mass of the exchange boson. For the case where two 
hadrons of equal mass m collide, these equations read in their cms 

q « v^/2, (1) 
qo < myJl- 4m 2 /s . (2) 

Implicit in these equations is, of course, a specific Lorentz frame with respect to which 
the energies are measured. In a simple bremsstrahlung experiment, these limits are easily 
satisfied by selecting only photons or dileptons of low energy in the laboratory frame. 
In the complex multiparticle systems formed in the course of nucleus-nucleus collisions, 
however, there are many binary collisions, and their respective cms frames do not generally 
coincide either with one another or with the overall nucleus-nucleus cms frame. 



In such a situation, it is better to look at relativistically invariant quantities, such as the 
dilepton invariant mass M. Looking at invariant masses means that go is no longer fixed 
but must vary over its full kinematic range, which for our example of colliding equal-mass 
hadrons is given by 

s- Am 2 + M 2 

In Figure 1, we show the three functions ([!]) and (j2|) and (|3|) for the case m = = 140 
MeV, my = m a ~ 500 MeV and dilepton invariant masses M = 10 and 300 MeV. It is 
immediately clear that the assumptions underlying the SPA are not fulfilled even for small 
M: the kinematic range accessible to go is never much smaller than the limits set by the 
SPA. The situation becomes even worse for larger M. 



3 Cross sections: approximate and exact 



In order to quantify what effect the use of the SPA has on yields, it is necessary to compare 
approximate cross sections to an exact formulation. All such dilepton cross sections can 
be written Q as the product of a purely leptonic prefactor 



K = {a/3ir) l + (2fi 2 /M 2 ) v /l-4 / u 2 /M 2 

(with // the lepton rest mass) times the cross section for production of a virtual photon 
7* with mass M, 

d(?hhi+£- = « da hhr . (4) 
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Now the traditional procedure has been to factorise the virtual photon cross section d<Jhh~i* 
into electromagnetic and strong pieces, writing it in terms of a current J M and three-particle 
phase space 



dR 3 = 5 4 (p a +p b -Px-P2- q)(d 3 p 1 /2E 1 )(d 3 p 2 /2E 2 )(d 3 q/2q ) 
dahh ^=4na^(-J»JJ\M h \ 2 ' Uh 



(5) 



where F is the incoming flux. It turns out that such factorization is unnecessary. A 
re-derivation yields exactly, without factorization j|, 

dM 2 



M 2 



dR 3 
(2tt) 5 F 



(6) 



where to, re run over all diagrams of the reaction irir — * tttt^*, including emission by the 
central blob. 




(GeV) 



Figure 1: The region between the dashed lines is the domain of integration for go when 
calculating cross sections as a function of dilepton invariant mass. The SPA is valid only 
when qo is much smaller than the two solid lines shown. This is clearly not the case. 

Since people do not like 3-particle phase space, dR% is usually simplified by neglecting q 
in the delta function, permitting reduction to 2-phase space dR 3 ~ dR 2 (d 3 q/2q ), with 
dR 2 = 5(p a +Pb-pi - p 2 )(d 3 p 1 /2E 1 )(d 3 p 2 /2E 2 ). Eqs. (@), © then simplify to the Riickl 
form 0, | 



j Ruckl 
uo hh£+e- 

dM 2 



k a 



2 M 2 



7T 



J da hh J\q^q 2 -M 2 J J%) , 



(7) 



where da^h is the on-shell cross section for the purely hadronic reaction irir — > irir. The 
current is first angle-averaged over the photon solid angle d£l q and then integrated over 
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the kinematic range M < qo < A = [s + M 2 — 4m 2 ]/2y / s of the photon energy </o- At this 
point an ad hoc factor to correct for the factorization of dR% is usually also included. 

A better derivation by Lichard || led to a form similar to Eq. (|7|) but without the (3/2) 
prefactor and with the inclusion of q and M in the current. 

In order to use the exact formulation (0), by contrast, it is necessary to formulate da^h^* 
in terms of three- particle phase space invariants, which for the schematic reaction a + b — > 
1 + 2 + 3 are defined [0 as s = (p a + p b ) 2 , h = (pi-p a ) 2 , si = (pi +P2) 2 , s 2 = (P2 +P3) 2 , 
and t2 = (j>b — Pz) 2 - The final-state phase space integral is then given by Q 

,„ , s 7T f dti ds 2 dsi dt 2 

dR 3 (s) 



4A 1 / 2 (s, to 2 , m 2 ) 

weighted by the 6x6 Cayley determinant B and where X(x,y,z) = (x 
One then obtains exactly 

k(M 2 )tt 



V 



j —exact 

aa hhe+e- 



dM 2 (2vr) 5 M 2 8\(s,m 2 ,ml) 



dt\ ds2 dsi dt2 



B 



52(M m nM* n 



nJtJ. 



(8) 
Ayz. 

(9) 



where all terms {M m )^{M* n )^ are written in terms of the five invariants. 



As a by-product of the 3-particle phase space language, one can define an intermediate 
approximation which, while still factorising the electromagnetic part out of da hh i+£-, 
writes the current J in terms of its invariants HI: 



da 



hM+e- 



4na 



kit 



dM 2 



(2tt) 5 M 2 8A(s, m 2 a ,m 2 ) 



dh \Mh(s,h)\' 



ds2 dsi dt2 



-J 2 (s,ti,s 2 ,si,t 2 ) 



(10) 



4 One boson exchange model 

Equation (|9|) may be exact in terms of the matrix elements A4, but it does not specify 
what A4 should be. For a quantitative comparison of the approximations (0) and (|l0|) 
to the exact cross section @, we must therefore turn to a simple microscopic model of 
pion-pion interactions Q. We use a simple one boson exchange model with the <r, p and 
/(1270) mesons included. The hadronic part of the lagrangian 

C = g a ad^-d^Tt + g p p^-irxd^ + g f f^ir-d v TT (11) 

is fitted to elastic tt + it~ collision data to fix the constants, then supplemented through 
minimal substitution by the corresponding electromagnetic interaction lagrangian. We 
work at tree level only. Monopole strong form factors are included for t and u-channels; 
no electromagnetic form factors are needed at present. Virtual photon emission is im- 
plemented gauge-invariantly for pion and p emission as well as the -rnra, irnf and imp 
vertices. 



This model is probably far from a perfect description of the pion-pion interaction, so that 
results obtained below can serve merely as a good pointer towards answering the question: 
How different are the results when we use the SPA or the exact cross section? 
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5 Numerical results 



To quantify the differences between the approximations and the exact formulation, we 
have studied the five distinct pion-pion reactions. Writing (+— ) — > (H — ) as shorthand 
for the reaction (tt + it~ —* ir + ir~£ + £~) and so on, we have calculated within our OBE 
model cross sections for (+-) -> (+-), (++) -> (++), (+-) -> (00), (00) -> (+-) and 
(+0) — * (+0). Numerical results were checked by performing a number of consistency 
checks, including testing for gauge invariance in the a, p and / sectors separately. 

Figures 2-4 show the cross sections da hh £+£- for dielectron production from the five distinct 
reactions as functions of ijs, for M = 10 MeV and 300 MeV respectively. Final-state 
symmetrization factors were included where appropriate. Initial-state symmetrization 
was also included for (+- 1-) — ► (++) and (00) — ► (H — ) in order to facilitate use within 
a thermal pion gas environment. Because they are identical in structure to their charge- 
conjugate versions, cross sections for the reactions (+0) — > (+0) and (++) — » (++) were 
doubled. All cross sections were computed using the same OBE model and parameter 
values. 



We see that there is a complete hodgepodge of curves, with a few underestimating the 
"exact" cross sections (solid lines) but most overestimating the exact curves by factors 
1-5. The largest discrepancy between approximations and the exact result occur for 
the reaction (++) — > (++) (Figure 4): for the Riickl approximation, factors 3 (for 10 
MeV) to 40 (for 300 MeV) arise, while the Lichard approximation yields corresponding 
overestimation factors of 1.9 and 14-20. 



In Figure 5, we attempt to cast some light on the physical origin of the discrepancies 
observed. Plotted are the cross sections for reactions (++) — ► (++) and (H — ) — ► (00). 
The upper lines correspond, as before, to the Riickl, Lichard and 3-phase space current 
approximations, while the solid line again represents the exact result. The lowest dash- 
dotted line, on the other hand, represents the exact result but excluding all internal 
diagrams and their cross terms with external ones. The difference between this lower 
line and the exact result (solid line) therefore represents the contribution of the internal 
diagrams; while the difference between the lower dash-dotted line and the upper lines 
(approximations) represents the change in cross section due to inclusion/exclusion of q in 
the externaZ-emission diagrams. We see that the contribution of internal emission is not 
all that large, albeit nonnegligible. By far the most important effect on da/dM is the 
neglect of the dependence of M on the photon momentum q. We believe that this neglect 
is at the heart of the considerable differences between approximations and exact results. 

Finally, Figure 6 shows the dilepton production rates per unit spacetime summed over all 
seven reactions, calculated using the Boltzmann formula p 



1 ds\{s,m ,m )——— (12) 



d 4 x dM 2 32vr 4 J v ' ' ' WsjT) dM 2 

for temperatures T = 100 and 200 MeV respectively. Right-hand panels show the corre- 
sponding ratios, obtained by dividing a given approximate by the "exact" result. Again, 

2 "Internal" diagrams are those for which the photon is emitted at a hadronic vertex or by the exchange 
boson. 
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the Riickl approximation is the worst, as expected; overestimation factors range from 2 to 
4 for T = 100 MeV, and 2-8 for T = 200 MeV. The Lichard and 3-phase space current 
approximations overestimate by factors 1.4-4, depending on temperature and M. 

Note that none of the approximations approaches the exact result for small values of 
M: even for the smallest value shown (M = 10 MeV), the discrepancy is still above 
40% for the Lichard and current approximations and larger than a factor 2 for the Riickl 
approximation . 



> 

CD 

o 



E 



b 



(+-)->(+-) 



!! M=10MeV 








0.004 



0.002 - 



(+-)-»(+-) 
M = 300 MeV 



s 1/2 (GeV) 








(GeV) 



> 

CD 



E 



b 



(+0)-»(+0) 

(-O)-k-O) 

M = 10 MeV 








0.004 



0.002 - 



(+0)-»(+0) 
(-O)-X-O) 

M = 300 MeV 



(GeV) 








(GeV) 



Figure 2: Cross sections for n + iv~ — > iv + iv~e + e~ and 7r + 7r° — ■> ■K + -K°e + eT for fixed invariant 
masses M = 10 and 300 MeV. Solid line: exact OBE calculation Eq. (|9|). Dash-dotted line: 
Riickl anproximation (pp. Dashed line: Lichard approximation. Dotted line: 3-phase space 
current Old). 



6 



-0.5 
> 

CD 

\ 
n 

E 



b 



(+-)-> (0 0) 
M= lOMeV 








0.002 



0.001 



(+-)-> (0 0) 
M = 300 MeV 



(GeV) 








(GeV) 



-0.5 
> 

CD 

\ 
n 

E 



b 



(0 0) ->(+-) 
M = 10 MeV 




0.002 



0.001 



(0 0) ->(+-) 
M = 300 MeV 








1 2 
s 1/2 (GeV) 







1 2 
s 1/2 (GeV) 



Figure 3: Same as Figure 2, for the reactions 7r + 7r — ► 7r°7r°e + e and 7r°7r° — > 7r + 7r e + e 
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Figure 4: Same as Figure 2, for the reaction 7r + 7r 4 
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Figure 5: Contribution of external-emission vs. external-plus-internal diagrams for the re- 
actions (H — ) -» (00) and (++) -» (++) for M = 300 MeV. Lines are as in Figs. 3 and 4. 
The new dash-dotted line below the (solid line) exact calculation represents contributions 
arising solely from emission of 7* by an external pion line, but taking q into account in M, 
in contrast to the approximations (upper lines) which neglected q. 
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Figure 6: Left: Total bremsstrahlung yields of e e for all seven pion-pion reactions, for a 
Boltzmann gas with temperatures T = 100 MeV and T = 200 MeV. Right: Ratios of SPA 
approximation calculations divided by exact OBE rate. Solid line: exact OBE calculation, 
Dash-dotted line: Riickl approximation, Dashed line: Lichard approximation, Dotted line: 
3-phase space current, including (upper) and excluding (lower) the mass of the virtual 7* in 
the current. 
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We hence believe that the SPA is flawed when used in a heavy ion context except in 
special situations. The suppression of the "exact" rates compared to traditional calcu- 
lations shown here would imply that bremsstrahlung dielectrons cannot make up for the 
discrepancies between measured dielectron rates and the cocktail of reactions used for 
comparison. Greater attention will have to be paid to the tttt annihilation, r/ Dalitz and 
other channels contributing at low M. The suppression we find would also have a bearing 
on calculations of the Landau-Pomeranchuk effect [pT[|. 
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